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We study the violations of Bell inequality for thermal states of qubits in a multi-qubit Heisenberg 
model as a function of temperature and external magnetic fields. Unlike the behaviors of the 
entanglement the violation can not be obtained by increasing the temperature or the magnetic field. 
The threshold temperatures of the violation are found be less than that of the entanglement. We 
also consider a realistic cavity-QED model which is a special case of the mutli-qubit Heisenberg 
model. 
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Historically the violation of Bell inequalities [l| was 
considered as a means of determining whether there is 
entanglement between two qubits. In 1989 Werner ||] 
demonstrated that that there exists states which are en- 
tangled but do not violate any Bell-type inequality, i.e., 
not all entangled states violate a Bell inequality. Further 
studies showed that the maximal violation of a Bell 
inequality does not behave monotonously under local op- 
eration and classical communications. So Bell violations 
only give a hint for entanglement and can be considered 
an entanglement witness Most recently Scarani and 
Gisin built the relations between Bell inequalities and the 
usefulness for quantum key distribution |p| and quantum 
secret sharing , and Diir showed that even multipartite 
bound entangled states can violate Bell inequalities ||] . 

An interesting and natural type of entanglement, 
thermal entanglement, was introduced in magnetic sys- 
tems and studied within the Heisenberg isotropic 
anisotropic Heisenberg models |10 1^ and the Ising 
model in a magnetic field |l2). Another natural type of 
entanglement was recently introduced in Bose-Einstein 
condensates The Gibbs state of a system at ther- 

modynamic equilibrium is represented by the density op- 
erator 



p{T) = exp {-H/kT) /Z, 



(1) 



where Z =tr[exp (— if/fcT)] is the partition function, H 
the system Hamiltonian, k is Boltzmann's constant which 
we henceforth take equal to 1, and T — 1/ f} the temper- 
ature. As p{T) represents a thermal state, the entan- 
glement in the state is called thermal entanglement 
At T = 0, p{T) represents the ground state which is pure 
for the non-degenerate case and mixed for the degenerate 
case. The ground state can be entangled but the thermal 
state p{T) dX T — oo can not be entangled as p{T) is a 
completely random mixture. The entanglement in the 
ground state of isotropic Heisenberg models was studied 
in the literature 15|. A recent interesting work fl^ ] 
showed that the success of the density matrix renormal- 
ization group method in understanding quantum phase 
transitions p7{ is due to the way it preserves quantum 
entanglement under renormalization. 



A complication in the analysis of thermal entanglement 
is that, although standard statistical physics is charac- 
terized by the partition function, determined by only the 
eigenvalues of the system, thermal entanglement proper- 
ties require in addition knowledge of the eigenstates. In 
other words thermal entanglement is determined by both 
eigenstates and eigenvalues of the density operator p{T). 
An interesting behavior of the thermal entanglement is 
that it can increase with the increase of temperature and 
magnetic fields P,p^. A natural question is if we can 
induce the violation of Bell inequality by changing tem- 
perature and magnetic fields. 

The study of thermal entanglement in the magnetic 
systems build a bridge between the quantum information 
theory and condensed matter physics. As there is inti- 
mate relations between the entanglement and the viola- 
tion of Bell inequality, it will be interesting to investigate 
the Bell inequality in the magnetic systems, specifically 
in some well-known solvable quantum spin models such 
as Heisenberg models. 

Arnesen et al. pointed out that one can test the 



Clauser-Horne-Shimony-Holt (CHSH) inequality |8|jl9| 
for the thermal state by measuring different components 
of the magnetic susceptibility tensor since the different 
components of the magnetic susceptibility tensor are pro- 
portional to spin-spin correlations in different pairs of di- 
rections. Very recently Mancini and Bose proposed 
a scheme for producing robust thermal entanglement be- 
tween atoms in distant cavities connected by an optical 
fiber, and the scheme serves as an experiment to detect 
thermal entanglement when the cavities are very near. 

In this paper we will study the violation of Bell in- 
equality for thermal states of qubits which interact via a 
multi-qubit Heisenberg interaction model. One particu- 
lar case of this model was recently realized in cavity-QED 
by Zheng |2|]. We will also discuss the violations in this 
realistic model. 

The most commonly discussed Bell-inequality is the 
CHSH inequality. The CHSH operator (a, a', 6, b' are unit 
vectors) reads 



B ^ a- a ® {b + b') ■ cf + a' ■ a ® [b -h') ■ a. 



(2) 
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In the above notation, the Bell inequality reads 



B) = tT{pB) < 2, 



(3) 



where p is an arbitrary two-qubit state. The maximal 
amount of Bell violation of a state p is given by |22] 



B = 2Vu + u, 



(4) 



where u and u are the two largest eigenvalues of TpT^, 
Tp being the 3x3 matrix whose elements is 



{Tp)nm = tr(/9f7„ (g) CT™). 



(5) 



Here ai = ax-,<J2 = o'y, and cts = Oz are the usual Pauli 
matrices. 

Now we define a quantity 



M = M + M - 1 



(6) 



which ranges from -1 to 1. The quantity M, which we 
call violation measure of Bell inequality in this paper, 
indicates the Bell violation when M > and maximal 
Bell violation when AI — 1. 

For any pure and some special mixed two-qubit state, 
there exists a simple relation between the concur- 
rence C (the measure of entanglement) and the max- 
imal violation of the Bell inequality: 



(7) 



We see from the above equation that the entanglement 
implies the violation of Bell inequality. But for a general 
mixed state, there is no this simple relation. 
We consider a 7V-qubit Heisenberg model 



N 



B 



N 



i=l 



(8) 



where J is the exchange coupling contant, A is the 
anisotropic parameter, and B is the magnetic field along 
z direction. In this model all particles interact equally 
with each other. The model reduces to the two-qubit 
model for N — 2 and to the three-qubit Heisenberg 
model with the periodic boundary condition for = 3, 
respectively. 

By using the collective spin operators Sa = 
(q^ — x,y,z), the Hamiltonian is rewritten 

as 



H = j{Sl + Sl + AS-^) + BSz 
= JS^ + J(A - l)Sl + BSz 



(9) 



up to a trivial constant. Due to the key term one can 
generate multipartite maximally entangled states in this 
system [pipll 



The partition function is obtained as ||ll|] 

N/2 N-2k 

Z = ^Nk ^-llJliA-l)im-N/2+kf]^-l3B{m-N/2+k) 
k=0 m=0 

'-■'N/2-k){N/2-k+l)^ ^^Q^ 



xe 



where Ni. 



0. 



N 

The reduced density matrix for any two qubits is given 
by 



Pi2 = i;+|00)(00| + + w(|01)(01| + |10)(10|) 

+y(|01)(10| + |10)(01|) 

= - [1 + {V+ - V^){aiz + (72z) 

+2y{aix(J2x + criyO-2y) + (1 - 4.w)aiza2z] (H) 
with matrix elements 

N^~2N + 4(S'2) ±A{Sz){N-l) 



v± 



AN{N -1) 



y 



4:N{N~ 1) ' 

2{Sl + Si) - N 



(12) 



2A^(A^-1) ■ 
The relevant expectation values are given by pl 

N/2 N-2k 

{f{Sz))^Y.^^ E f{^-N/2 + k) 

k=0 m=0 

^^-rJ.Il{A-l){r,i-N/2+kf]^-PB(m-N/2+k) 



xe 



-f3J{N/2-k){N/2-k+l) 



/z, 



N/2 N-2k 

{Si + si) - E ^'^ E - k){Nl2 -k + l) 

k=0 m=0 

-{m- N/2 + k)'^] 



^g-/3.7[(A-l)(m-JV/2+fc)ig-/3S(m-W/2+fc) 



xe 



-l3J(N/2-k)(N/2~k+l) 



/z. 



(13) 



Here f{Sz) is an arbitrary function of Sz- Eqs.(|lOD-(|l^ 
gives completely the reduced density matrix. Along the 
standard procedures for calculating the maximal viola- 
tion and the concurrence, they are obtained as 

M EE M{ J,A,(3,B) 

= 8y^ + (1 - 4u;)2 - min[V, (1 - 4^)^] - 1 (14) 
C = C(J,A,/3,i3) =2max{0,|2/|-VT^} (15) 

For TV = 2, Eqs.(|l|) and (|l|) become 

M = {2sinh2(/3J) + [cosh{(3B)e-'^^-^ - cosh(/3J)]2 

- min{sinh2(/3J), [cosh{l3B)e-^^^ - cosh(/3J)]2}} 
/[cosh{(3B)e-l^^-^ + cosh{(3J)]^ - 1, (16) 



C = max{0,sinh(/3|J|) 



-0AJ 



} 



/[cosh{l3B)e-'''^'^ +cosh{f3J)]. 



(17) 
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From the analytical expressions of M and C we imme- 
diately obtain 

M(J,A,/3,B) = A/(-J,-A,/3,B), 
C(J,A,/3,S) = C(-J,-A,/3,S), (18) 

M(J,A,/3,i?)-M(J,A,/3,-B), 
C(J,A,/3,S) -C(J,A,/?,-i?). (19) 

From Eq.(p^ we see that the violation of BcU inequal- 
ity is invariant under the simultaneous sign changes of 
the parameters A and J. The entanglement also has the 
same property. In particular the violation does not de- 
pend on the sign of J when A = which implies that 
the violations are same for the antiferromagnetic (J > 0) 
and ferromagnetic cases (J < 0) in the XX model. From 
Eq.(|l9|) we know that M are symmetric with respect to 
B. Although this conclusion is obtained for the case of 
N — 2, it is valid for any N in our model. 

Fig.l gives numerical calculations of the violation mea- 
sure and concurrence as a function of temperature for 
different magnetic fields. The concurrence is plotted in 
this paper for a comparison. It is clear from Fig.l that 
there exist threshold temperatures Tm (Tc) for the vio- 
lation measure (the entanglement). After the threshold 
temperature Tm (Tc) the Bell inequality is not violated 
(the entanglement disappears). Note that Tc is inde- 
pendent on the magnetic fields, while T^/ decreases with 
the increase of magnetic fields. We also observe that 
Tm < Tc which implies the fact that some entangled 
states do not violate a Bell inequality. In the temper- 
ature range Tm < T < Tc, the state is entangled but 
the Bell inequality is not violated. Now we see the case 
of B = 2.5. The entanglement can be increased by in- 
creasing the temperature, however the Bell inequality is 
not violated at any temperatures. The comparison shows 
that the behavior of the violation measure is very differ- 
ent from that of the concurrence. 



FIG. 1. The violation measure and the concurrence as a 
function of temperature for different magnetic fields. Viola- 
tion measure: -6 = (cross points), B = 1 (circle points), 
B — 2 (diamond points) and B = 2.5 (box points); The 
concurrence: B — (solid line) , B = 2 (dashed line), and 
B — 2.5 (dotted line). The parameters J = A = 1 and N = 2. 

In Fig. 2 we plot the violation measure as function of 
temperature for different number of qubits. The results 
show that although we can have entanglement the Bell 
inequality can not be violated for TV > 3 as seen clearly 
from the figure that M < 0. When iV > 3 we obtain 
the two-qubit reduced density matrix by tracing out all 
the other qubits but two of them. This usually increases 
the entropy and make the Bell inequality much less be 
violated. We also see that the threshold temperature 
Tc increases as the number of qubits N increases and 
Tm < Tc- 

Then we see if the magnetic fields can induce the vio- 
lation for > 3. The violation measure versus the mag- 
netic field is given in Fig. 3. The temperature is choose 
to be close to the absolute zero temperature. From Fig. 3 
we see that we still can not obtain the violation of Bell 
inequality for A'^ > 3 by increasing the magnetic fields 
although the entanglement can be increased by the mag- 
netic fields. For the case of A = 2 we find the threshold 
magnetic field Bm for the violation measure and Be for 
the concurrence. The threshold magnetic field Bm is 
less than Be- We also find that the threshold magnetic 
field increases as the number of qubits increases. The 
magnetic field can increase the entanglement. The more 
complicated behaviors of the entanglement can be found 
for the case of 6 qubits (dotted hue) . There are two dips 
when the magnetic field increases. After the two dips, 
the entanglement increases and then decreases and reach 
the threshold point. 
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FIG. 2. The violation measure and the concurrence as a 
function of temperature for different number of qubits. The 
violation measure: N = 2 (cross points), A'' = 3 (circle 
points), and N = 5 (diamond points); The concurrence: 
N = 2 (solid line) , iV = 3 (dashed line), and TV = 5 (dotted 
line). The parameters J = — 1 and B — A = 0. 



From Eq.(|ll|), the violation measure D is obtained as 

D = v+ log2 v+ + log2 w_ + (w - y) log2(w - y) 
+ {w + y) log2(w + y) - (v+ + w) log2(u+ + w) 
- (u_ + w) log2 {v-+w). (23) 
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FIG. 3. The violation measure and the concurrence as a 
function of the magnetic field for different number of qubits. 
The violation measure: N — 2 (cross points), A*' = 3 (cir- 
cle points), and A'' = 6 (diamond points); The concurrence: 
iV = 2 (solid line) , iV = 3 (dashed line), and iV = 6 (dotted 
line). The parameters J — —1, T = 0.05, and A = 0. 

Another interesting feature of quantum entanglement 
is that a bipartite entangled state (systems A and B) 
can be more disordered locally than globally j26| . If we 
measure the disorder by von Neumann entropy S{p) = 
— tr(p log2 p) , one inequality holds for all separable states 

S{A),S{B)<S{A,B). (20) 

which is called disorder inequality in this paper. For the 
models we are considering, S{A) = S{B) holds. For the 
general case S{A) is usually not equal to S{B). Then we 
define a quantity 

D = S{A)~S{A,B), (21) 

which can be larger than zero and acts as a quantitative 
measure of the violation of the disorder inequality (|20| ) 
Interestingly the quantity D gives directly a lower bound 
of the entanglement of formation Ef{A,B) of the state 
|3,i.e., 

Ef{A,B)>D, (22) 
which shows that the state is entangled ii D > 0. 



FIG. 4. The violation measure D of the disorder inequality 
, the concurrence C and violation measure M as a function of 
T. We plotted C (solid line), M (dashed line), and R (dotted 
line) for N — 2 and C (cross points), M (circle points), and 
R (diamond points) for A*' = 3, respectively. The parameters 
J = -1 and B = A = 0. 



The numerical results are given in Fig. 4. We find that 
the behaviors of the violation measures are similar. For 
N = 2, both the disorder inequality and Bell inequality 
can be violated. But the threshold temperature Tu of 
the violation measure is different from Tm and Tc, i.e., 
Td < Ti\i < Tc- There exists a temperature range Td < 
T < T]\j in which the Bell inequality is violated while 
the disorder inequality is not violated. For A^ = 3 we 
still can not find any violation of both the disorder and 
Bell inequality although the state is entangled for a range 
of temperature. 

Finally we consider a realistic model in cavity QED. 
Consider the interaction of A^ qubits with a single ca vity 
mode in the vacuum state via Tavis-Cumming model psj . 
For the case of large detuning 6 between qubits and the 
field mode, the effective Hamiltonian is given by 

He^X{Sl + Sl + S,), (24) 

where A — g^/S and g is the atom-cavity coupling 
strength. Comparing Eqs.(p^) and we find A = 
and J = B = X. Note that A can be negative due to the 
detuning. The numerical results are shown in Fig. 5. We 
see that the entanglement can be obtained while the Bell 
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and disorder inequality can not be violated. The behav- 
ior of this physical model is similar to that obtained from 
the 'abstract' model. Both the analytical and numerical 
results show that the Bell inequality is more difficult to 
be obtained than the entanglment. 
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FIG. 5. The violation measures and the concurrence as a 
function of A/T. We plotted C (solid line), M (dashed line), 
and R (dotted line) for N — 2 and C (cross points), M (circle 
points), and R (diamond points) for A'^ = 3, respectively. 



We consider the multi-qubit thermal state in the multi- 
qubit Heisenberg model which is exactly solvable. We 
trace out all the other qubits but two. Then we ob- 
tain the two-qubit reduced density matrix, from which 
we have studied the violation of Bell inequality. The re- 
sults show that the Bell inequality are violated for the 
case of = 2 but not violated for the case of iV > 3 at 
any temperature and magnetic fields in our model. Un- 
like the behaviors of the entanglement, the violation can 
not be induced by increasing temperature and magnetic 
fields. In comparison with the thermal entanglement, 
they are relatively hard to be obtained. 
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